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Abstract. In this work, we introduce a new discretization to the frac- 
tional Laplacian and use it to elaborate an approximation scheme for 
fractional heat equations perturbed by a multiplicative cylindrical white 
noise. In particular, we estimate the rate of convergence. 
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1. Introduction 

In this work, we are interested in the space approximation of the solutions 
of fractional stochastic heat equations. Equations where leading operators 
are fractional or more general pseudo-differential operators are widely used 
to model complex phenomena. For example, they are ubiquitous in the 
study of the quasi geostrophic flow, the fast rotating fluids, the dynamic 
of the frontogenesiaj in meteorology, the diffusions in fractal or disordered 
medium, the pollution problems, the mathematical finance and the transport 
problems, see e.g. [21 El EJ dTJ [19], [26J ETJ [32] and the references therein. The 
wellposedness of these equations, in the deterministic and stochastic cases, 
has been extensively studied see e.g. [21 [H EJ EJ El [TBI ED]- Although the 
numerical approximation of the solution is needed in applications, the num- 
ber of numerical schemes relevant to such approximations is quite restricted. 
The main difficulty of the numerical approximation of fractional equations 
is related to the fractional operator. For example, contrarily to second or- 
der differential operators, the fractional operators can not be discretized by 
three points. Using the classical schemes and as a global operator, all the 
values on the grid should be used in every step. Via the integro-differential 
representation of the fractional operator, a direct discretization is based on 
the discretization of the integrals. This idea has been used in the numerical 
study of the deterministic conservation law driven by fractional power of 
the Laplacian in [10J. Unfortunately, as it is mentioned in the paper, the 
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convergence of the scheme elaborated is slow and leads to some unreason- 
able values. Moreover, in the theoretical study, the author did not give an 
explicit form to the discretized operator. The discretization of the integrals 
has been already used for the Liouville and Riemann fractional operators 
and has yielded the Griinwald formula, see for short list [JBJ G3] • 

During the preparation of this work, we found the work of Westphal [31J, 
on an approach to fractional powers of operators via fractional differences. 
The approximation given by Westphal provides a rigorous mathematical 
support to a numerical discretization of fractional operators. Westphal de- 
fined the fractional operator of the infinitesimal generator A of a semigroup 
S(t), defined on Banach space X, by: 

A r := s - lim T~ r (I - S(t)Y, < r < 1, 

where 

(I - S(r)Y := £ C-^S(jr) = £ — ^ S(jr) 

and s — lim means the strong limit, i.e. the limit in the space X of t~ t (I — 
S(t)) v /, for every f £ X. In particular, she proved: 

[x/t] 

(1.1) lim T~ r C J T ' X fi x ~ i T ) = DT f( x ) 

3=0 

where [x/t] is the integer part of x/t and D r is the fractional differential 
Riemann-Liouville operator defined on R + , An intuitive way to discretize 
the operator D r can be obtained by taking r = ~: 

[nx] 

To encounter the difficulties of the direct discretization of the fractional 
operator, probabilistic technics have been used. In particular, in [25] the 
authors used the Monte Carlo method to approximate numerically the so- 
lution of some deterministic fractional partial differential equations, among 
them the Burgers equations. In [22J the authors used wavelet techniques to 
approximate the Kolmogorov equation driven by the infinitesimal generator 
of a Feller process. 

Our idea to discretize a fractional operator A" 1 is to discretize first the 
operator A then to take the fractional power of the discrete operator. As 
far as the authors know, this idea is new. 

In this work, we discretize the fractional Laplacian, in the way described 
above, and we elaborate a scheme to approximate the fractional stochastic 
heat equation. Our aim is also to calculate explicitly the rate of the conver- 
gence and to show its dependence on the fractional power of the Laplacian. 
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We are also interested in the critical values of the fractional order which 
insures the convergence of the scheme as well. 

We consider the following fractional stochastic heat equation: 

( fu(t,x) =£zu(t,x)+g(u(t,x))^(t,x),t>0,xe(0,l), 

(1.2) I u(0,x) =u (x), xe (0,1), 
[ u(t,0) = u(t,l) = 0, t > 0, 

where -J^ = (—A) a", a > 1 is the fractional power of the Laplacian. Let 
us denote by A a = Aa, where A = —A is the Laplacian with boundary 
Dirichlet conditions, defined on D(A) = H 2 ' 2 (0, 1) n H^' 2 (0, 1). H k 'P(0, 1), 
for k G N, p G [1, oo) is the Sobolev space of order k. The fractional operator 
A a = A? is well defined, see e.g. Lemma 2.6.6 in [20J and it is given by the 
formula (see [20] pp 72-73): 

(1.3) A a = (-A)f := / z%- l A(Iz + Ay l dz. 

K Jo 

The operator A a is a closed densely defined operator with domain of defini- 
tion given via the complex interpolation of order D(A a ) = [H, D(A)] s , 
see e.g. \15\ I29j and Theorem 4.2 in |28| . More precisely, 

oo 

D(A a ) = D{A a ' 2 ) = {v€ L 2 (0, 1) : £ Agug < oo}, 

fe=i 

where = (v, e^) = \/2 w(x) sin/c7rx<ix and A& = k 2 ir 2 , G N are the 
eigenvalues of the operator A corresponding to the eigenfunctions: {e^ = 
\/2 sin A;7r-}fc g N. The map g : R — >• R, is a bounded Lipschitz continuous 
function on R. The operator g is regarded as a nonlinear operator from 
H = L 2 (0, 1) to C(H), the set of bounded linear operators on H, defined by 
g{u){h) = {(0,1) 9x1—7- g{u(x))h{x) G R}. In other words, the nonlinear 
operator g is the Nemytski map associated with function g. For v G H, 
g{v) is given as a multiplicative operator. From the hypothesis that g is 
bounded we have ||<7(i>)|| < &o ; where &o = sup K |g(x)|. {W(t),t > 0} is a 
cylindrical Wiener process on the probability space (fi, J 7 , {J}t>o, IP). The 
initial condition uo is a L 2 (0, l)-valued J-o-measurable function. In section 
we will suppose stronger condition on the diffusion term g and on the 
initial condition uq. In particular we will suppose that A s g is a bounded 
Nemytski map for some 5 > 0; ||^4 <5 5(v)|| < b$ with 6^ = sup K lA'^x)! and 
uq belongs to a given fractional Sobolev space. 

We rewrite the equation (jl.2p in the following form: 

, , r d-u(i) = -A a u(t) dt + #(«(*)) dW(t), t > 0, 

( ' I «(0) =n . 
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Let us denote by {S a (t),t > 0} the semigroup generated by — A a and by 
H v ' 2 the fractional Sobolev space of order v. By a solution of the equation 
(11, 4j) , we mean, see e.g. [8]: 

Definition 1.1. Suppose that a > 1. An J-t-adapted H 8,2 -valued continuous 
process u = {u(t),t > 0} is called a mild solution of equation ftl.4\ ) with 
initial condition uq G H 71,2 ,!] > 0, iff for some p > a 2 { t _ 2 e where 9 < 
min{^yi, rf\ 

(1.5) E sup \u(t)\ p g2 < oo, T>0 

te[o,T] 

and for all t > 0, a.s. the following identity holds 

(1.6) u(t) = S a (t)u + f* S a (t - s)g(u(s)) dW(s). 

Jo 

We introduce the space: 

Definition 1.2. Let T > andp G [1, oo] be fixed and H a Hilbert space. By 
Zt,6,p(H) we denote the space of all H e,2 -valued continuous and J-t-adapted 
processes u = {u(t),t G [0, T]} such that 

(1.7) IMIr,fl,p := E SU P = IMIiPff^fo.TJxife. 2 ) < °°- 

te[o,T] 

If 9 = 0, we use shortly the notation Zt p (H). 

The following result of existence and uniqueness of the solution of equation 
(|1.4p can be concluded from the calculus in [1]: 

Theorem 1.3. Let a > 1 and let uq be a H v ' 2 -valued J-Q-measurable func- 
tion such that 

^\uo\ P Hv , 2 < oo 

for some p > ^ry, < 9 < min{^i — ^,rj} and let T > 0. Then there 
exists a unique mild solution u G ^T,e,p, of equation \1.4\) - 

The paper is organized in the following way. In section [2] we describe the 
discretization of the fractional operator. In particular, we apply the idea for 
the Galerkin approximation and for the finite difference method. In section 
Owe elaborate a numerical scheme to approximate the solution of the frac- 
tional stochastic heat equation (II. 4p . In section HI we give some preliminary 
estimations of the Green functions corresponding to the fractional operator 
and to the approximated operator. The section [5] is devoted to prove the 
convergence of the approximated solution to the solution of the equation 
(|1.4j) . In the end of this introduction, let us mention the following refer- 
ences where the the approximations of certain stochastic partial differential 
equations are treated [U [T2l HH HU [21] . Let us also mention that, we take in 
the whole paper p, a > 1 and the values of the constants may change from 
line to line. 
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2. Discretization of the fractional operator 

Let us first recall notions about the approximations of the Laplacian —A, 
see e.g. p]. We consider the Gelfand triple (V, H, V'), where V H = 
H' c — >■ V densely, where V = Hq(0, 1) and V is its dual. The operator A 
defines a coercive bilinear form on V by 

a(u,v) := (Au,v) = (u,v), u, v E V, 

where u', v' are the first derivatives of u and v in the distribution sense. It 
is widely accepted that to approximate the Laplacian —A, it is sufficient to 
approximate the bilinear form a. 

Let V n be a finite dimensional subspace of V generated by a basis (/j)j =1 . 
It is easily seen from the formula: a(u n , v n ) = Y^j=i u i v j a (fii fj)-, f° r u v„ '■= 
Y^ r j=i u jfj-> v 'Vn, := Y^ r j=i v jfj e that the projection of the bilinear form 
a on V ra , denoted by ay n := aP\; n is well defined via the double series index 
(a(fi, fj))i,j- Using Riesz representation we can rewrite the bilinear form 
ay„ as: 

n 

«V„ (Uv n , VV n ) = (A Vn U Vn , V Vn ) = ^ U i V i ( A V n fi ,fj), 

where Ay n is a positive bounded linear operator on V n which is well defined 
via the stiffness matrix A Vn := (a(fi,fj) = {A Vn f iy fj))i >j=1 . 

Now we define the fractional power of the approximated operator Ay n by 
the following formula (see [20] pp 72-73): 

- sin7r^ f°° a , , i 
(2.1) A$ n = j o z~ l A Vn {zI + AvX 1 dz 

and the fractional bilinear form 

a 

a a. 

The fractional stiffness matrix A-^ is then given by (ay njQ ,(/j, fj) = (Ay^fi, fj))fj = 
Our idea is to investigate how and on what rate the operator A^ and the 

bilinear form ay n)Q! are good approximations to the operator A a respectively 

the bilinear form a a := (Azu, v). 

Before going through this calculus, let us apply this method to calculate 

the approximation of the fractional operator and of the stiffness matrix 

corresponding to the following two methods; the Galerkin method and the 

finite difference method. 

2.1. Approximation by Galerkin method. Let V n be the subspace gen- 
erated by the basis (ej)" =1 defined above. Recall that ej(.) := \^2sm(jir.). 
It is known that the operator A and Az are diagonizable under the basis 
(ej)JL 1 . Consequently the approximating operator A\> n and the stiffness 
matrix Ay n are diagonal with respect to the basis (ej)" =1 and with corre- 
sponding eigenvalues (Aj)f =1 . Thanks to the boundness of the approximation 
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operator Ay n , it is also easy to define Ay n on the basis as: Ay ei = A?ej, 
see also Lemma IC.ll 

2.2. Approximation by finite difference method. Let {xi := ~,i = 

0, . . . , n} be the set of grid points and let (pi, i G {1, . . . , n — 1} be a pyramid 
function, i.e. a function, which takes value 1 at the grid point -, vanishes 
at the other grid points and is linear between the grid points. The approx- 
imating space V n _i is then generated by {(^i^Zi ■ For an implementation 
reasons, we will focus on the stiffness matrix. It is well known that the 
stiffness matrix Ay n ^ corresponding to the finite difference approximation 
of the operator A is given by: 

'2, i = j, 

-1, i = j + 1 or i = j - 1, 
0, otherwise. 

Using the matrix theory it is easy to calculate the fractional power of the 

a 

Av„_i, denoted by Ay via the formula: 
sin §7r / _ , 



(2.2) A^ = — <— j z^Av^ilz + Ay^yUz. 



3. Discretization of the fractional stochastic Heat equations 

Let us first observe that the eigenvalues of Ay n _ x are given by Xj n '■= 
j 2 ir 2 c jn , j = 1,2,. ..n - 1 where c jn := sin 2 (|^)/(^) 2 , the corresponding 
eigenvectors e™ = ((ejk)k) are given by pQ: 

ejk = \ -sm(j -7t). 
V n n 

From Lemma lC.lt it is easy to see that A ? j = 1, 2, ...n — 1 are the eigenval- 
ues of Ay corresponding to the eigenvectors e" = ((w ^ sin(j^7r))fc), j = 
1, 2, ...n — 1. The semi group S™(t) generated by Ay is given by: 

n-l 

(3.1) S n a {t)x := {(Y,Gl(t,i,j) Xj ),), 

a 

where G™(t,i,j) := X]fc=i e ~ Ufcne fcj e fcj and x = (%j)i<j<n-i- Let us de- 
fine now the operators: P n : L 2 (0, 1) -> M n_1 and E n : R™" 1 -)• L 2 (0, 1), 
called projection respectively interpolation operators given by the following 
formula: 

For all / G L 2 (0, 1) and for all x G R n_1 : 

n-l 

(3-2) P n f := J2( p nf) k e%, 

k=l 
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where 

n— 1 

(3.3) (i^)k e^fo). 

i=i 

rt— 1 

(3.4) £ n z:=^<x,e£)e fc , 

fe=i 

It is easy to see that the operators P n and E n satisfy the properties: 
Lemma 3.1. • P n and E n are bounded linear operators such that 

1 1 -Frail < 1 an d \ \En\ \ = I- 

• P n E n = I n , where I n is the identity matrix in E™ -1 . 

• E n P n = P n , where P n is the projection on the finite dimensional 
space in I? generated by ei, e 2 , e n _L i.e. P n f :=Y^Zl{f,ej)ej. 

• P n ej = e™ , if j = 1, 2, ...n — 1 and zero z/ j > n. 

• £ n e™ = e^, j = 1,2, ...n - 1. 

• {ej, j = 1,2, ...n— 1} are i/ie eigenf unctions of the operator E n Ay P n 

a 

corresponding to the eigenvalues X? n . 

• The Green function of E n Ay P n > which is also the kernel of the 

semigroup: E n e Vn ~ 1 P n , is given by 

n-l a 

C%(t,x,y) :=Y,e- tX ^e k (x)e k (y). 
k=l 

Lemma 3.2. For < 6 < \ + fa and 7 > - — 4-, £/iere exists K > 0, suc/i 

— 4 4 ' a a 7 " 

#wrf V/ £ (0, T], 

ct a 

\S a (t) - E n e- tA Vn-ip n \ c(H ^ D{A _ S)) + P~ 5 (5 Q (t) - ^ne^-iPn)!!^ 

(3.5) < K<f) at s,j{t,n), 
where 

(3.6) 4> aAl {t,n) :=( 



fjr n 1 + a— 4<5 , -1 

n id +n at 2a , i < <5 < A + | Q . 



Let us make the following convention to write 4> a ,s^ shortly as 4>a,S when 
7 is not presented, i.e. when \ < 5 <\ + fa. 

Let us now discretize the diffusion term g. We denote by g n the matrix 
which is given by the column vectors: (P n ((goE n )ej))j , 1 < j < n — 1, where 
o is the composition of the two operators: E n and the Nymetsky map g. For 
y 6 E n_1 , the operator (go £7 n )y acts as the Nemytski map associated with 
the function g{E n y{.)), i.e. ((g o E n )y)h = {(0, 1) B x i-> g(E n y(x))h(x) G 
E}. We denote by W„(t) the vector (B^t), B 2 (t), £ n -i(*)) of independent 
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Brownian motions. We introduce the following multidimensional stochastic 
differential equation, where u n (t) := (Wn(*))l<A<n-l : 

( 3 7 \ J du n (t) = -A$ n i u n (t)dt + g n (u n (t)) dW n (t), t > 0, 

\ u n (0) := P n (u ). 

For the existence and uniqueness of the solution of the stochastic differential 
equation (|3.7p . we refer e.g. to Theorem 2.3 in [TT] : 

Theorem 3.3. There exists a continuous-W 1-1 valued Ft — adapted process 
Un = {u n (t),t > 0} solution of the problem {3. 7| ) such that : 

r-t 



(3.8) u n (t) = e ^-in„(0)+ / e { * s)A ^-i g n (u n (s))dW n (s), a.s. 

Jo 

Furthermore, there exists a constant CV,n, such that 
(3.9) 

E sup |u n (i)|^„_ 1 < C Tn ||g||(l+E|?xo|^ 2 ), for all T>0 and pe[l,oo). 
te[o,T] 

We define the I?— valued stochastic process u n (t) := E n u n (t). We prove, 
see Appendix [S] that: 

Lemma 3.4. The process u n {t) := E n u n (t) satisfies the following stochastic 
integral equation: 

u n (t) = E n e- tA ^-iP nUQ + / E n e- (t ~ s)A ^P n g(u n (s))dW n (s), a.s, 

Jo 

(3.10) 

where 



ft _ § n-l „ t 

/ E n e- {t - s)Av ^P n g(u n (s))dW n ( S ) := V / E n (e 
Jo j=1 Jo 

= V / E n e- [t ~ s)A ^P n g(u n (s))e 3 dB,{s) 
3=1 Jo 

and (e * s ^ Av n-i g n (u n (s)))j is the j 's column of the matrix e s - )A v„_i g n (u n {s)) . 
Furthermore, for all T > and p £ [1, oo), 

E sup \u n {t)\l 2 <C T ^ m {l + nu \ P L2 ). 
ie[0,T] 

In other words, u n (t) satisfies the stochastic partial differential equation: 
f du n {t) = -S„A| i i P n n n (t) dt + g(u n {t))dW n {t), t > 0, 

j u n (o) =E n p nUo n . 1 
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4. Preliminary estimates 

In this section we give a priori estimations to the Green functions, G™ and 
G a , corresponding to -^ — A? respectively — Ay and to their difference. 

Lemma 4.1. For < 5 < \ + |a and for all 7 > ^ — 4^, i/iere exzsis 
K > 0, suc/t fftotVt G (0,T], 

n—l 

(4.1) ^A- 25 |e-' A , f -e-'4| 2 < i^n~ a t _1_ « +4 ^ 



00 



(4.2) ^A- 2i e- 2 ", Sif^-^t^^W+n-"^!,,^,^ 

where Xb is the characteristic function of the set B: 

1,6 e 5 

Proof. To get the estimation (|4.ip . we use the mean value theorem. We 
obtain, 

n — l a a 71—1 a 

Xj 2& \e~^ 1 - e~ X ^ l \ 2 = ^ Xj 25 e- 2t{wj)a \l - e^'^^ 1 | 2 
j=i 3 1 

< ^A7 2 V M W)°|l-c| l | 2 (j7r) to * 2 e ^ 1 -4)^) a ^ )T G [0,1] 

n— 1 a 

3=1 

We know that: < c| n := | sin(g)/(g )\ a < 1. Taking 1 - c| n = 

0(g) 2a , we obtain 

n—l a a n—l 

^AT 25 |e- A /'-e- A -f < Kt 2 n- 4a J2f a ~ 45 e- 2a+1 i at 
3=1 3=1 



n-l 



< ^ 2 n- Q ^/ Q - 4<5 e- 2Q+1 ^* 

/»oo 

< ^ 2 n- a / x 3a - 4<5 e- 2Q+1 ^*dx 





/■oo 

< ^- 1 -i+ 4 ^n" a / y 3a - 45 e- 2Q+1 ^dy. 
J 
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The integral f^°y 3a 4(5 e 2a+1 v a dy converges provided < 5 < \ + fa. 
Hence, there exists K > such that: 



n— 1 

(4.3) 

To get the second estimation (|4,2p . let us first consider the case 5 6 [0, |]. 
We use the known result: for all 7 > 0, there exists a constant k = k^ such 
that e _a: < if ;pf, we get 



j=n j=n (2A ? t) 7 A 2<5 j =n 



But ET=n(j) a7+i5 < K Ji + °° x- a ^- i5 dx < 00, provided ± -4| < 7. Hence 



(4.5) J2 X i 25e ~ 2Xft ^ Kn a ^h-\ 

j=n 



For 6e(H + fa), we use the inequality: Y%Ln XJ 2S e~~ 2X i 1 < E^AJ 2 *, 
than we arguing as above and using the condition 8 > j ? we obtain 



OO q: OO OO 

(4.6) X>T 25 e- 2A 7* < ^AT 25 <Kn- 45 ^(-) 4<5 <^n- 45 . 

j"=n J=n i=n ^ 



□ 

Lemma 4.2. For < <5 < I + fa and 7 > — — 4—, there exists K > 0, suc/i 

— 4 4 'da 7 7 

i/iaf Vi G (0, T]j we /iave 



(4.7) sup \A- s (G a {t,x,.)-G n a {t,x,.))\ H <K4> aAl {t,n), 

xe[o,i] 



where <j> aj s^{Jb,n) is given by $3.6}) . 
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Proof. Using the definitions of the functions G a (t,x,y) and G™(t,x,y) and 
the fact that the orthonormal basis (ek)k>i C L°°(0, 1), we get 

\A- s (G a (t,x,.) - G n a (t,x,.))\ H 

Tl 1 a a OO a. 

= I E( e " A/ * " e- X ^)A- s e 3 (x)e 3 + £ e~ X T * A~ & e^e^ 

j=l j=n 

TL 1 a a OO a, 

< I E( e " A/ * " e-^AjefaMH + \ ^2e~ x ? t A- s e j (x)e j \ H 

j=l j=n 

< I E( e ~ A7 ' " e~ X ^ t )Xj 5 e,(x)e J \ H + \ e~ X l "xf e 3 (x)e 3 \ H 

3=1 3=n 
i=l i=n 

72 1 ct a OO a 

i=l i=n 

(4.8) 

Replacing P~2j) and (f¥7T|) in we get the result. 

□ 

As a consequence of Lemma 14.21 we obtain: 
Corollary 4.3. Under the same conditions in Lemma \4--S\ 

(4.9) \A- s (G a (t, ., .) - G£(t, ., .))|hxh < ^ ai< 5 j7 (t,n). 

O: 

Proof of Lemma 13.21 Let f € H. The semigroups 5 a and E n e Vn ~ 1 P n 
are acting on an element / via their Green functions in the following: 

(S a (t)f)(x)= I G a (t,x,y)f(y)dy 
Jo 

and 

(E n e- tA l-iP n f)(x) = J GZ(t,x,y)f(y)dy. 
Applying the Holder inequality, we get: 



o 
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A- s (S a (t)-E n e~ tA ^P n )f\ H = [ \A- 5 ((S a (t)-E n e' tA ^P n )f)(x)\ 2 dx]^ 

Jo 

= [[ \[ A~ s (G a (t,x,y)-G n a (t,x,y))f(y)dy\ 2 dx]^ 
Jo Jo 
1 r l 



< [ dx( \A- 5 (G Q (t,x,y)-G n a (t,x,y))\ 2 dy)]>\f\ H 

Jo Jo 

< [[ \A- s (G a (t,x,.)-G n a (t,x,.))\ 2 H dx] 1 *\f\ H 

Jo 



< \A- s (G a (t,.,.)-G n a (t,.,.)\ HxH \f\ H . 

Using (|4.9p . we get the result. For the second estimation, we have by a direct 
application of the definitions of Hilbert-Schmidt norm and the properties of 

a 

the semigroups S a (t) and E n e Vn - 1 P n , 

a CO Ol 

\\A- s (S a (t)-E n e~ tA ^-iP n )f HS := ]T \A~ 8 (S a (t) - E n e~ tA ^ P n ) ej \ 2 H 

Tl 1 a a OO a 

i=l 3=n 

Tl 1 ex, a_ OO a 

j=l j=n 

(4.10) 

Thanks to the basic inequality (a + 6)2 < a 2 + 62 for a, 6 > and by the 
formula (14. D and (I4.2p . we get the result. 



5. Convergence of the scheme 
Now we are ready to give the main result of this work. 
Theorem 5.1. For a > 1, 

la la 
4 4 4 4 

, . ,11 a, -13 

(5.2) max {__ + _} <£< - + -a 

and 

2a a 2a 

(0.0 w > maxj , — , — ^ 

v ; y 1 a-2'25-l'8<5-a-2 / 

assume that: 
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• (H g ) the diffusion term is the Nemytski operator defined by the Lip- 
schitz function g G D(A S ) and such that bs := sup xgR |A*5>(a;)| < oo. 

• (H UQ ) the initial condition uq is an D(A n )— valued LP random vari- 
able i.e. n G L p ^l,D{A^)). 

Then u n = {u n (t),t G [0, T]} converges to u := {u(t),t G [0, T]} in the space 
Zt, p (L 2 (0, 1)). Furthermore, there exists a constant K > 0, such that, 

\Ht)-u n (t)\\ ZTAL 2 M < K nuQ \ D{An)M n~^ 

(5.4) 

where £ is given by 

(5.5) e = min{|, 25}. 



In particular, for 1 < a < 2, the rate of convergence £ :- 



2 ' 



Remark 1. Let us remark that the rate of convergence is independent of 
the regularity of the diffusion term when the dissipation order a is less than 
the Laplacian dissipation. In this case, it is enough to take g G H 1 . 

Theorem 5.2. Assume that a > 2, p > g and uq satisfy respectively 
(Hg) with 5 = and (H UQ ), with n satisfying (15. ip . Then u n = {u n (t),t G 
[0, T]} converges to u := {u(t),t G [0, T]} in the space Zt :P (L 2 (0,1)) and 

\Ht)-u n (t)\\z T , p (LHo,i)) < KTMowtorr^-*-®. 

(5.6) 

First let us introduce some lemmata which we will use in the proof of the 
convergence. 

a 

r -a" 2 " t 

Lemma 5.3. The operator A commutes with S a (t) and with E n e Vn ~ 1 P n , 

for all t > 0. 

Proof. For the proof see the Appendix B. □ 

Lemma 5.4. Suppose that z G L q (0, 1) for q G [2, oo]. Let Z denote the 
multiplication operator by z. Then, for (3 < 1 — — , there exists a constant 
K > 0, such that 



/ s-e\\E n e- sA ^P n Zf HS ds < K(Y>-^)) 
Jo k=i 



z\\ q < CO. 



(5.7) 



Proof. Let us first estimate the term \\E n e sAy ™-i P n Z\\ 2 HS . Using Lemma 
13.11 we have 
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oo 

\ "2" ^ < — oA^ ^ « — ^ _ o A "2" 



\\E n e sA ^-i Pn Z\\ 2 HS = J2\E n e sAv r,-, Pn Ze k \ 2 L2 =J2\ZE n e 8Av n-x Pn e k \l 2 

k=l k=l 
n—1 a n—1 a 

= Y, \ZE n e- sA ^e n k \ 2 L2 = Y \ze- x ^E n e n k \l 2 
k=i k 

n-1 a 

= /Z e ~ X * nS \ ze k\h- 



Let us observe that by the Holder inequality, | | ^2 < where 
\ + ^ = \- Moreover, since \&k\h 2 = 1 an d |efc|L°° = 2 1 / 2 it follows by 

applying the Holder inequality that \ek\h r < 1 x l q . Let us recall that A^ n := 

(cfc„(/c7r) 2 )f , where c kn := ^J^ , we get, 

V In. I 



n-1 



k=l 

n-1 

-(|c fcn |fc7r) Q s 



< 2«|z|| 9 ^e 

fc=l 

Therefore and thanks to the fact that: (f ) a < | Cfcn |§ = | sin(^)/(^)| a < 1 

/* OO tv f* OO ^ " " 

Jo Jo fc=1 

= 2i|z|£, V(c fcn )2(^- a )(fc7r)^- a / i 
fc=i ^ 

-foo poo 

< K\z\\ q V(fcvr) a/3 - a / T- p e^drr 
k=i Jo 

+oo 



fc=l 



Since by our assumptions a(l — /3) > 1 the series on the RHS above is 
convergent and the result follows. □ 

The following Lemma is a special case of Lemma 2.7 from [4j: 

Lemma 5.5. Provided that v > p v the operator R u : L p (0, T;L 2 (0, 1)) — > 
C([0,T];L 2 (0,1))) given by 

R u h(t)= [ (t- S y- l S a {t- s)h(s)ds, h G L p (0, T; L 2 (0, 1)) 
Jo 
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is well defined, linear and bounded. Moreover, there exists a constant C PtU > 
such that for all h G L p (0, T; L 2 (0, 1)) 

(5.8) \ R uh\ C ([0,T];L 2 (0,l)) - K p,uT U P Hlp(0,T;L 2 (0,1))- 

Lemma 5.6. Let j < 5 < \ + |ct and p > max{l, a ^L^g }, then there 
exists v , satisfying max{0, ^ + 2a — 77 + |}<^<l? such that the operator 
G n>v : LP(0,T;D(A 5 )) -> C([0, T];L 2 )) given by 

G n>v h(t)= / (t-sy- 1 [S a (t-s)-E n e~ it ~ s)Av n-i Pn ]h(s)ds, h G L p (0, T; D(A 5 )) 
Jo 

is well defined, linear and bounded. Moreover, there exists a constant K p> t > 
such that for all h G 17(0, T; D(A 5 )) 

(5.9) \G n ,vh\c([o,T};L2) < K T , v {n' 2h + n~f ) | h\ LP ^ T . D{AS)) . 

Proof. Let us fix h G L p (0,T; D(A S )). Then for t G (0,T) and thanks to 
Lemma 15.31 we have 

\G n , v h{t)\ L 2 = | [\t- S y- l [S a (t-s)-E n e- {t - s)A ^-iP n ]A- 5 A s h(s)ds\ L 2 
Jo 

< / (t- s r- 1 \A- 5 [S a (t-s)-E n e- {t - s)A ^-ip n ]A 5 h(s)\ L 2ds 
Jo 

< / (*-*r- 1 [|5 a (t- a )-^ n e- (t -' ) V-ip ft [| |A*fc( a )| La 



o 

From Lemma l3T2l ( (|3.5p and (|3.6p ) and by applying the Holder inequality, 
we get 

\G n ^h{t)\ L 2 < K l\t - «)"-V«,i(* " «, n)|A*fc(a)| £a 



o 



< 



T „ p _i rT 



J 



Thanks to the basic inequality: (x + y) e < cq(x 6 ' +y ), for 6* > 1 and x, y > 0, 
we have 



l' T {s u - l (t> a , 5 {s,n))T^ ds = [ T (s u -\n- 2S + n- i ls~ 1± ^ M ))^ ds 
Jo Jo 



Up 

(5.10) 



o r p _ a p / / -J 1 + a — 4(5 \ p 
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The last integral in the RHS of (|5.10p converges provided v > ^ + 2^ — 77 + ^ • 
Hence 



T 



( / (/" 1 ^(s,n)) p!l ^) Ppl <c p , T (n~ 2<5 +n-t) 





The choice of v such that f + 77 + |<^<li s possible thanks to the 

2a 



condition p > ^df+iir Finally, we have for all t > 

\G u h(t)\ L 2 < K T ^ p (n- 25 + n'^)\h\ Lv ^ T . D(AS)) . 

Proof of Theorem 15.11 Let 

(5.11) M n := \\u - u n \\ p T = E sup - u n (i)|^. 

te[o,T] 

We have from equations (|1 .6|) and (|3.10p . 

Ht)-u n (t)\ p H <c p (A(t)+B(t)) 



where 



(5.12) 



A(t) := \S a (t)u -E n e tAv ^P n u \ p H 



B(t) := I f S a (t-s)g(u(s))dW( S )- f E n e (t s)K ^P n g(u n {s))dW n (s) 
Jo Jo 

(5.13) 

Estimation of A{t): Using Lemmata [53] and E21 f (j33|) and We ob- 

tain, 

a 

A{t) := \(S a (t)-E n e' tA ^P n )u \ p H 

= |A-"(5«(t) - Ke'^-Puji^ol^ 

< |(S Q (i) - J E n e-* A ^ 1 p n |^ ( ^ D(A _ t))) | J 4^o|^ 
< 



l n °l_D( J 4 ? ') ( ^a,'7^' n ) 



Thanks to the condition: rj > I + j, the power — 1+ ° Ar] is positive, hence 



A(t) < K T |u |^Jn- 2 " + n-f) p 
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Consequently: 

E[supA(t)l < K T E\u \ p D(Av Jn-^ + nSy 

[0,T] v ; 

< K T E\u \ p D(Av) (n- 2 ^ + n~ S r). 

(5.14) 



Estimation of B(t): Let us first introduce the transformations defined 
on the set of Nemytsky maps N, such that for h £ N: 

n 

(5-15) h(u)e r .= { 'i^y^ 

Then we write the second stochastic integral in RHS of (|5.13p . as 

/ E n e~ {t ~ s)Kv ^P n g(u n {s))dW n {s) := V / Ene"^-^ Pn g( u n {s))ejdBj 
Jo j=1 Jo 

t a 11 

= ( E n e- {t - s)K ^P n ^?Jsr)dW{s). 
Jo 

(5.16) 

Using the factorization method see e.g. [1] and [8], we can again rewrite the 
integrals in f)5. 13|) for < v < 1 as: 

t i-t 

v-l 



S a (t- s)g(u(s))dW(s) = / (t- sY^Sait- s)Y{s)ds 
o Jo 

and 

a n a 

f E n e- {t - s)k ^P n 'glJ?{sf)dW{s) = f\t - sy^E^^-i P n Y n (s)ds, 
Jo Jo 

where 

(5.17) Y(s):= [ (s-ry u S a (s-r)g(u(r))dW(r) 

Jo 

and 

n 

(5.18) Y n (s) := f\s - r)-»E n e- A ^ s -^P n 'gls?l^dW(r) 



o 
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Consequently, 



B(t) 



{t-s)' / - 1 S a (t-s)Y{s)ds- / (t-s) u - L E n e 



-In -(*- 



-ip n y n ( s )ds|^ 



< c p (\ / (t- S y- 1 S a {t-s)(Y(s)-Y n (s))ds\ 



''-I- L ' / / .1 // „ (' S ) A V„_ 1 p 



+ I / (i-s)"" 1 S a {t-s)-E n e 
Jo 1 

(5.19) 

Using Lemmata 15.51 and 15.61 and taking 



Y n (s)ds\ p H ) 



, _i 1 1 2<5 1. 

maxjp , - + h -} < ^ < 1 

2 za a p 



we deduce that 



E 



[0,71 



[0,T] JO 



su Pj B(t)j < c p Esup(| / (t-s) i/ " 1 5 a (t-s)(y(s)-y n (s))ds|^ 



+ | / (i-s)"" 1 5 a (t-s)-£; ne 



y n (s)ds|^ 



< K T>P [E\Y - Y n \l P{QtT;L2{QA)) + (n- 2 * + n-t f E|y n r LP(0jT;D{ ^ 



(5.20) 



Calculation of E|y n |^ T . D r^s\y By the Burkholder's inequality and Lemma 
15.41 there exists a constant C p > such that 



T 



E\Y n (s)\ P n(AS ,ds<C p / E( / (s-r)- 2 »\\A 6 E n e [s r)A ^ P n g(u(r))\\ 2 HS dr 



o x Jo 



r 



< K p l\( (\ s -r)~ 2u \\E n e 



^P n A- 5 A S g(u{r))f HS drY ds 



JO 
T, r s 



^PuWhs su P \ a 9{u{r))\ Lo , {H) dr) ds 
re[o,T] ' 



<K p b p 5 



r-^\\E n e '^-iP n \\ 2 HS dr)Us 



10 JQ 

-oo 



< K p b p s T^k- a{1 - 2 " Y ) ' ■ 

k=l 

(5.21) 



Since p > max{ 2 £ i _ l , ^rj} and 5 > ^, we can chose v, such that 

, _i 1 1 1 25, 11 

maxj p , — | 1 f < z/ < 

lF ' p 2 2a a ; 2 2a 

what implies that a(l — 2z^) > 1 we infer that the last term in (I5.21j) is finite. 
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Calculation of E|F—Y n |^ 0T . L 2( i))- Using the formula (|5.17j) and (15.180 . 

we have: 

\Y(s)-Y n (s)\ p H < c p (| f\s-T)^{S a {s-r)-E n e- {s - r)k ^P n )g{u{r))dW{r 

n 

P n {g(u(r)) -^J?)))dW(r 



+ 



(s - r)- v E n e (s r)A ^«-i 



By Burkholder's inequality, we have 

E\Y(s) - Y n (s)\ p H < c p (e(J\s - r)- 2 »\\{S a (s - r) - E n e~ {a ~^»-i P n )g(u(r))f HS dr 



+ E(J ( s -r)-^\\E n e 
(5.22) 



2;-n ;.- ,~( s - r )A^ i _ 1 



P n {g(u(r))-g(u n {r)))\\ 2 HS dr 



Using the well known functional inequality ||A.B||#s < H^H^sl-Bl^^, we 
estimate the first term in the RHS of (|5.22j) as follow: 



\s r)- 2u \\{S a {s-r)-E n e {s r)k ^P n )g{u(r))\\ 2 HS dr 



< / (s-r)- 2 »\\A- d (S a (s-r)-E nt 
Jo 



<S r)Av ^P n )A s g(u(r))\\ 2 HS dr 



< I [ s - r ) 
/o 



- 2v \\A- 5 [S a {s-r)-E n e 



^P n )\\ 2 HS \A s g(u(r))\ 2 dr 



< b] ( s -r)~ 2 »\\A- 6 {S a (s-r)-E n e (s r)A ^ P n )\\ 2 HS dr 
Jo 

(5.23) 

Thanks to ([375]) ([375]), we have 



(s r )~ 2u \\A- s (S a ( S -r)-E n 



- (S - r)A ^P n )\\lsdr 



< K / {s-r)- 2v 4> 2 aS {s-r,n)dr 

Jo 

< K(n- 45 I r- 2v dr + n- a I' r'^-^+^dr 



) 



The integral f£ r 2v dr is finite thanks to the condition v < 



and the integral L s 

25 1_ 1 

a 2a' 2 2a 



2 2a ^ 2 

s r _2l/_1_ ~ + ~(ir converges thanks to the conditions v < 



min{— — 7T-, i — tt-}. Hence, we take the parameter ^ which satisfies the 
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following inequalities: 

(5.24)max{p \ - + - + - - - } < u < mm{- - - - -}. 

The parameter v exists thanks to the conditions: 5 > | + |- and p > 
max{ ggjfb^ , 2^1 }■ Hence 

(* - r)- 2 "[|5 a ( S -r)-£? n e- ( '" r)A? »-iP n [|^ s dr 
< K T (n- 4<5 + n- a 



o 



(5.25) 

By replacing (I5.25j) in (I5.23p . we get 

a 

(s - r)- 2 "||(5 a ( S -r)-£; n e- ( - r)A ^- 1 p n ) y ( u ( r ))||2 rs d r 



o 



Hence, 

- r)- 2 "||(5 a ( S -r)-£7 ne " ( '" r)A? »-iP n )i7(u(r))||? rs d'-)' 

(5.26) 



Now we estimate the second term in (|5.22p . Arguing as in the proof of 
Lemma 13.21 we have 



\E n e {s r)A Vn-i P n (s(u(r))- || 



2 

HS 1 



n 

^ | ((/(«(r)) - 9 (u») )^e- (s - r)Av - 1 P^. 
i=i 



From Lemma l3.lt we have 

^^P.ej = I e- (s - r)A -e,, j G {1, • • •, n — 1} 

1 0, j > n, 
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n 

using the definition of g(u n (r)), the Lipschitz property of g and Lemma[3J] 
we get 

n 

a 



||^e ^ r ^»-i Pn (sMr))- 5 (u>)))||^ 

n-l a 

= 53 e" 2( "" r)A £ I (s(u(r)) - s(«"(r))) ej\ 2 H 



3=1 

n—1 c, 

< Y, e ~ 2iS ~ r)Xl \9(< r ))-9(u n (r))\l {H) 

3=1 

n—1 a 

< ^2e- 2{s - r)X ^\u(r)-u n (r)\ 2 H . 
3=1 

n-l a 

< sup \u(r) -u n (r)\ 2 H V e - 2(s_r)A ^. 
r6[0,a] J = 1 



Hence, 

n 

e(£(s - rr^ll^e-^^-iPn^CtiCr))-^^)!!^) 1 ) 



n-l 



P 



< E sup \u{r)-u n (r)\ p H ( V / r~ 2iy e- 2rA /»dr) 

re[Q,«] j=l 

(5.27) 

Arguing as in the proof of Lemma 15,41 we get a constant K which depends 
only on u, such that 

n 

OO £ 

< ifE sup mr)-u n (r)|^fVr a(1-2v) V 

re[o,»] y 

< sup |n(r) -u n {r)\ p H 
re[0,s] 

(5.28) 

since < f < | — Now replacing (15.26h and (15. 28ft in (j5.22j) . we obtain 

P s IT _. f f„-4S i „-a\ 2 , it? „„„ „,n/ MP 



E|F( S )-r n ( S )|^ < 2^,^ n-^ + n-^r +E sup \u(r) - u n (r)\ p H 

VV 7 re[0,s] 

(5.29) 
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Fl~K — V \ p 

^l 2 lrl lLP(0,T;L 2 (0,l)) 







E\Y{s)-Y n (s)\ p H ds 



' K„ r /„ ( ( >r 4S + n~ a ) 2 + [ E sup |u(r) - u n (r)|^ds 

' h re[o,s] 



< K p ^ h Jn- 2S P 



+ n 2"+ / E sup \u(r) — u n (r)\ p H ds 

JO re[0,s] 



(5.30) 



and by replacing (|5,2ip and (|5.30[) in ()5.20f) . we obtain 



E 



sup B{t) 

L [0,T] 



p,T,b 



n~ 25p + + I E sup |u(r) - u n (r)|^ds 

re[0,s] 



(5.31) 



Finally, we join the estimations in (|5.3ip and (|5,14p . we obtain, 
E[ S u V \u(t)-u"(t)f H ] < K,t,\u \ d(AV) J [ T M SU P l«(r) 

[0,T] ' UID ( A ')' V Jo re[Q,s] 



(5.32) 

Thanks to Gronwall Lemma, 



E[sup|u(t) -u n {t)\U < K p Tl . 

L j ^' w wli/J — p,T,\u \ D(AV) , 



Now it is easy to see that from the conditions on rj, that % < 2rj, hence 



E[sup|u(i)-^(t)|^] < K P nuolD(Avhb n-Z p , 



(5.33) 
where 

26, which implies that £ = §■ 



where £ is given by (|5.5|) . Furthermore, if 1 < a < 2, then ^<^ + ^<l< 



□ 



To prove Theorem 15.21 we will use the same scheme as before, but with 
different estimations: 
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Lemma 5.7. For 7 > < j3 < 1, there exists K > 0, such that Vi G 

(o,n 

\S a (t) - E n e- tA ^-ip n \ c{H) + \\S a (t) - E n e~ tA ^-iP n \\ HS 

(5.34) < K(n- a 2t-2+n ~t 35—), 

Proof. The proof follows the same steps as in Lemmata 13.21 and 14.11 and by 
considering the following key estimation: 

n— 1 a a n—1 

£ |e-V* - e- A 7"*| 2 < Kt 2 n- 4Q ^i 6Q e- 2Q+1 ^ 
j=i i=i 

n-l 

i=i 

POD 

< Kt 2 n- a+al3 / x 3Q - Q/3 e- 2a+1 ^*^ 



Jo 

(5.35) 

□ 

Lemma 5.8. Provided p > ~ < 7 < and | + ~ — 1 < /3, i/iere 

exists v satisfying: max{ 1+ ^~ a ^ + ^> 2 + pi ^ ^ ^ 1> suc ^ that the operator 
G n , u : LP(0,T;L 2 (0,1)) -> C([0,T];L 2 )) given by 

G n , v h{t) = / (t-s)^ 1 ^^-^-^"^"^-^^^^)^, ft£ L p (0,T;L 2 (0,l)) 

is well defined, linear and bounded. Moreover, there exists a constant K PtU > 
such that for all h G L p (0, T; L 2 (0, 1)) 



(5.36) |G n)!/ /t| ca0iT] . i2) < K p ^(n~ a 2 + n ~)\h\ LP{0)T . L 2^y 
Proof. Let h G L p (0, T; L 2 (0, 1)). Then for t G (0,T) we have 

\G n>u h(t)\ L * < / (i-sr-l[5 a (t- S )-E n e" ( *" s)Av n-ip n ] 11^1^)1^2 ds. 







From Lemma 15.71 and by applying the Holder inequality, we get 



1 ol — o.0 1 + a — cm/3 N 







\G n!U h(t)\ L 2 < K I (t- s)"- 1 (n~ a 2 (t - s)a +n a (t-s) —)\h(s)\ L2 ds 

f*T f*T 



, f 1 11 a-aB 1+a-ctfl . P . P- 1 . f 

< K( (s u - l (n- a 2 S -l +n 2^s —))—^ds)~{ \h{s)\ P L , 

Jo Jo 
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Thanks to the basic inequality: (x + y) 6 < co(x e +y e ), for 9 > 1 and x, y > 0, 
we have 



1 _ 7 _7 a — g/3 1 + a — a/3 V n , 7 P / /,, -i 7 ^ V 

(s v (n "5s 2+ n —2— s 2^ ))p-i ds < c p (n "a P -i / 2>P-ids 



q — a/3 



+ n 2 p-i 1 

JO 

(5.37) 

For p > 2 ^" 1 ; ~ < 7 < 2^- there exists ^ > 2 + p' sucri that the first 
integral in the RHS of (|5.37p converges. Furthermore, if/3 > | + ^ — 1 
then we can chose v satisfying also the inequality v > + | and 

hence the second integral in the RHS of (|5.37p converges also. Hence for 
v > max{ 1+ 2~ al3 + |, | + |}, we have: 



/ -, _1+o-oJ>j_ 
V 1 ^ - 1 - 2a ' p-1 



p-1 



, I -I I I U-t _L — p Lt Lt f-r f ^ £ p <_£ t_t [J 

(/ (s (™~ Q 2 s ~2 + n r~ s 2S p <c p (n~ a 2+n r~ ) 



Hence for all t > 



a — a/3 , 



|G„/l(t)| £ 2 < ^T,p(n Q 2+n 2 )|/l| L p (0 ,T;L2(o,l)) 



□ 



Proof of Theorem 15. 21 We arguing as in the proof of Theorem 15.11 We 
define 



(5.38) 

M n := \\u-u n \\ p T =E sup \u(t)-u n (t)\ p H <c p (E sup A(t)+E sup B(t)), 

te[o,T] '' ^ te[o,T] te[o,T] ' ' 

where A(t) and B(t) are given respectively by (|5.12p and (|5. 13|) . The term 
Esup tG r 0i T] A(t) is estimated thanks to the inequality (|5.14p . To estimate 
the term B(t), we use Lemmata 15.51 and 15.81 

ft 

e' 



supS(i) < CpEsup I / (t-s) u - 1 S a {t-s)(Y(s)-Y n (s))ds\ p H 

L [0,T] J [Q,T] ^ Jo 



+ I / (t-s)"" 1 S a (t-s)-E n e 
Jo 



y n (s)ds|^ 



< K T,p 

(5.39) 



E|y-yj? 



CK — Ctj3 



n lLP(0,T;L2(0,l)) 



LP (0, T;L 2 (0,1)) 



Now we calculate E|Y - ^1^(0^^2(0,1)) and E I^Ilp(o,t ; l2(o.i))- By the 
Burkholder's inequality and Lemma [5 41 there exists a constant C p > such 
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that 



^E\Y n (s)\ P L2 ds < K p jE(J\ 8 -r)- 2v \\E n e (s r)K ^P n g{u(r))\\ 2 HS drY 2 ds 
< K p [ T E( /"(a-rr^H^^^Pnll^blioodr)^ 



JO 

(5.40) < K p b p Q J o [J\- 2v \\E n e' r ^P n f HS drf ds 

+ 00 p 

< KXT(Y, k ' a{1 ~ 2V) Y ■ 
fc=i 

Since < v < \ — Jj- what implies that a(l — 1v) > 1 we infer that the last 
term is finite. In the aim to get an estimation to E|Y — Y n \ P LP ^ Q T . L 2( ]\)> we 
use the inequality (]5.22p . Let us remark that the estimation (]5.28p remains 
true for the second term in the RHS of (|5.22p . Let us now estimate the first 
term in RHS of this inequality. We have: 



/< 

J o 



S s r )-^\\S a {s-r)-E n e {s r)K ^P n g{u(r))\\ 2 HS dT 



< I ( s -rr 2 »\\S a (s-r)-E n e P n \\U\9(<r))\\ 2 * 



„,-2ni c / ,.-(•— '-'m- , d i2 I,,., ,, .mi 

C(H) f 

< b 2 { s -r)- 2v \\S a (s-r)-E n e- {s - r) ^-,p n f Hsdr 
Jo 



(5.41) 

Thanks to Lemma 15.71 we have 



\s r)- 2v \\S a {s-r)-E n e {s r)k ^-iP n f HS dr 



o 



< K n 



-a+a/3 



The last two integrals converge provided v < min{4 — 4~, ^-o^}- Hence, 



E (j Q ( s r)- 2v \\{S a {s-r)-E n e {s T)A ^P n )g{u{r))\\ 2 HS dr)' 

■D 



(5.42) 
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By accumulating the conditions, the parameter v should satisfy: 

A+a-a8 1 7 1, . c 8 II-7I 1, 

maxj 1 — , — I — \ < v < mm] , , \ . 

1 2a p 2 p J l 2 2a 2 '2 2a J 

(5.43) 

Since a > 2 and p > it is possible to choose 7 £ - |) and 

a + | + (notice that p >-^^I<i-I^-L + i + I< 
1). With such choice of 7 and /?, we can find v satisfying (|5.43p . In fact 
^ + i- l<^ + ^<i + i + i</3<land7 + i + ^<i + i + i</3. 

a p a p a p 2 " ' a p a p 2 " 

From (pT28|) and (f5Ti2|) . we have 



supB(i)l < K TiP Jn- a ^ +a ^ + n~ a ^ + [ E sup \u(r) - u n (r)\ p H ds 
[0,T] J L Jo re[0,s] 

(5.44) 

Finally, we join the estimations in (|5.44p and (|5.14)) and applying the Gron- 
wall Lemma, we get 



\u-u n \\ T ,p := E[sup|u(t) -u n (t)\ p H ] 



1 



' [0,T] 



(5.45) 



To get a good precision of estimation, we push 8 to its lower bound and 7 
to its upper bound, i.e. we take 8 := i + -7 + \ and 7 = ^ — ^ then we get: 

\\u - u"\\ T , p < K nuo]DlAV)ibo + n-i + „"«f +«f + „-! <W 

< j{„ 1 , fc f n -2 W + n -(f4-f-f)\F 

- - fV r,|«o|D(A'7),bo\ v ' t T ^ 

(5.46) 

Thanks to the condition (|5.ip . we have ^ — § — | < | + ^ < 2?]p. Now, 
it is easy to get the inequality (|5.6|) . 

Appendix A. Proof of Lemma 13.41 



We apply on the both sides of equation (|3.8f) the operator E n , we get 
(A.l) 

E n u n (t) =E n (e- tA ^u n (0)) +£„( y e~ MA ^-i fl „K(s))^ n (s)). 

Using the definitions of the two operators P n and i? n and by the fact that 
P n E n = I, the RHS in (jA"7[) is equal to u n {t). Moreover 

a_ a 

E n e~ tA ^u n (0) = E n e~ tA ^P n u . 
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Let us now explain here how to get the stochastic term. We denote by (•)& 
the component of a vector and by (-)kj the component of a matrix. We have 

n-l nt 



E n (J^ e {t s)A ^-i g n (u n (s)) dW n {s)) :=Jl{f o e ^ 9n(u n (s))dW n { S )) ^ 

E/ Ev^'^- 1 ^^))), .dBiWefc 
fe=i ^° j=i 3 

n— 1 „in~l a 

e y E( e-(t_8)Avi ^K( s )))^^i( s )- 



j=i fe=i 



But 

n-l 



£(e ( ^ )Av "-^nK( S ))) fc / fe = ^n(e '"Vift,^^)))., 
fe=l 

where (e ^ s ^ Av ™- 1 (7 n (n n (s)))_ ; . is the column "j" of the matrix e ^ 8 Ay »-i g n (u n (s)). 
Hence we have the first result: 

ft y n — 1 a 



We know from the basic calculus on matrices that: 

By the definition of the matrix g n , the column "j" of the matrix g n (u n (s)): 
((9n(u n (s)))j) is equal to P n ((g o E n )ej), hence we have 



( e -(* s)A v.- l5n ( Un ( s ))) = e -<* s)A ^( 5n K( s ))), 



V / S n e- Av "-i (i - s) P ri <7(n"( S ))e^ J B,( S ) 
~T ■/ o 



3=1 

and we denote these integrals by 



/ ^e- Av "-l ( *" S) P n5 K( S ))^( S ) 
JO 



To get the estimation of u n (t), we use Lemma 13.11 and Theorem 13.31 then 
we have 

E sup \u n (t)\ P L2 <E sup H^nllM*)!^-! < C T>n iu||(l+EKI^ 2 ), for each T > 0. 

te[o,T] te[o,T] 
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Appendix B. Proof of Lemma 15.31 

a 

In fact, the operators S a (t), E n e Vn ~ 1 P n , and A are bounded and we 

a a 

have A- s S a (t)ej = S a (t)A- s ej and E n e~ A ^-^ P n A- s ej = A~ & E n e~ K ^-i t p n6j , 
for all ej,j G N. 

Appendix C. Lemma ICTT1 

Lemma C.l. Let (\k)k>o be the sequence of the eigenvalues corresponding 
to the eigenf unctions (e k )k>o of the positive operator A. Then {e k )k>o ar ^ 

at — 

also eigenf unctions of A 2 corresponding to the eigenvalues (\ k )k>o- 
Proof. Using the definition of the fractional operator (11.30 : 

A%e k = ^_ / t^~ 1 A(tI + A)' l e k dt 

v Jo 

■ Q'7T /'OO 

= / tv- 1 \ k (t + \ k )- 1 dte k 

n Jo 

(C.l) 

By the residues theory, we have 

(l-e im ) / ^- l {t + i)- 1 di = 2niRes{-l,^- l {t + £,)- 1 ) = -2irie™^ 
Jo 



so. 



/"OO 

JO sm ~T~ 



□ 
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